A Scaling Approach to Ideal Quantum Gases 
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The thermodynamic properties of ideal quantum gases are derived solely from dimensional argu- 
ments, the Pauli principle and thermodynamic relations, without resorting to statistical mechanics. 

PACS numbers: 



I. INTRODUCTION 

In this article we present a derivation of the ther- 
mal equations of state of ideal quantum gases by using 
solely dimensional analysis, thermodynamic considera- 
tions and simple physical arguments 1 !. The thermody- 
namic relations of gases depend strongly on the energy- 
momentum relation of the gas particles and on their spin, 
i.e., whether the particles are fermions or bosons. The 
energy-momentum relation of free massive particles is 
given by 



c 2 p 2 



mc 2 . 



(1.1) 



m denotes the rest mass of the particles, c is the velocity 
of light. Here we have subtracted the rest energy mc 2 
from the particle energy which results in a redefinition 
of the chemical potential, i.e. fj, has to be replaced by 
fi—mc 2 further on. In the non-relativistic limit mc 3> |p|, 

r> 2 

E p w whereas in the ultra- relativistic limit mc <C 
|p|, E p pa c|p|. Both cases can be written in the single 
exponent dispersion relation 



Z± ~Jv IPI 



i 1 -^ 2 -", (1.2) 



where v = 1,2 for relativistic and non-relativistic parti- 
cles, respectively. 

Another group of particles to which our considerations 
apply are quasi- particles in condensed matter systems. 
In these systems deviations from the state of perfect order 
can be described by a superposition of quantized elemen- 
tary excitations which behave as bosonic particles with 
an energy momentum relation also of the form 1)1. 2f) . The 
quasi-particles we are considering, as well as photons, 
have a vanishing chemical potential, i.e. their number is 
not fixed. At low temperatures this number is small and 
hence their interaction can be neglected. Quasi-particles 
are characterized by a fixed value of v > — . 7^ has to be 
expressed in terms of the quasi-particle parameters. Ex- 
amples for the case v — 1 are lattice vibration of solids 
with acoustic phonons as quasi-particles^ (c denotes now 
the sound velocity) or spin waves in anti-ferromagnets 
with magnons as quasi-particles^ (c denotes here the spin 



wave velocity). An example with v = 2 are magnons in 
ferromagnets^, where m ~ h 2 /(b 2 T c ). b denotes the lat- 
tice constant and T c the Curie temperature. For further 
examples of quasi-particles see e.g. reference*. 

It is interesting to remark that our considerations can 
be extended to non-integer values of v. One example 
are capillary wave excitations of the surface of superfluid 
i He, the so-called ripplonsi which show a dispersion re- 
lation l|1.2|) with m/c ~ ph/a and v = 3/2. Here a and 
p denote surface tension and the liquid density, respec- 
tively. Our results apply also to this case, but in the 
further discussion we will mainly focus on the cases of 
integer values of v. 



II. DIMENSIONAL ANALYSIS 
A. The classical case 

Considering a homogeneous gas of N point particles 
confined in a container of volume V = vN and coupled 
to a bath of temperature T, the free energy can be written 
in the form 



F(T, V, N) = NF(T, v, 1) = Nf(T, v) 



(2.3) 



f(T,v) is the free energy per particle. The mean sepa- 
ration a between particles is related to the volume per 



particle v 



where d denotes the dimension of the 



system. The pressure follows then from p — ~(df /8v)t- 
We begin with the consideration of a classical gas with 
fixed exponent v of the dispersion relation (|1.2|l . In this 
case the pressure can only depend on the three param- 
eters 7„, a and T. There is no further parameter. The 
temperature T appears in thermodynamics in the rela- 
tion T _1 = 8S/8E, where S and E denote the entropy 
and the energy, respectively. Entropy is determined only 
up to a multiplicative (Boltzmann) constant^. For the 
moment we will only assume that this constant is dimcn- 
sionless, i.e. we measure T in energy units. 

Denoting the dimension of mass, length and time by 
M, L and t, respectively, the dimension of 7^, a, T and p 
are M x ~ v L 2 - V t v - 2 , L, ML 2 t~ 2 and ML 2 ~ d t~ 2 , respec- 
tively. Since the three parameters have different dimen- 
sions, there is a unique combination of these which has 
the dimension of pressuraS. 
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as expected. We fix now our temperature scale by choos- 
ing the constant of proportionality in (|2.4|) to be equal 
one, i.e. pv = T, which is the known result. 

This suggest f(T,v) ~ Tin v. Clearly v should occur 
in the free energy in combination with another volume to 
make the argument of the logarithm dimensionless. This 
shows the limitations of classical physics. 

In the case of massive particles there is only a depen- 
dence on m or c as long as we use the single power 
law relation l|1.2|l with v = 1 or 2 fixed. However, 
for the description of the cross- over between the non- 
relativistic and the relativistic regime both m and c oc- 
cur and hence a new dimensionless parameter T/(mc 2 ) 
appears. The relation (|2.4|l is formally correct even in 
the region T ^> mc 2 . However, one has to expect that 
in this region pair creation processes will take place and 
a grand canonical description is more appropriate, as we 
will consider in section III CI 

In the case of charged particles the ideal gas properties 
are only preserved under certain conditions. If we con- 
sider for simplicity electrons (on a neutralizing positively 
charged background) the electron-electron interaction is 
of the order e 2 /a. To have still an ideal gas the Coulomb 
energy has to be small compared with the thermal energy 
i.e. e 2 /a «Tora T = e 2 /T < a-ii. 

B. Quantum case 

In the quantum case the Planck constant of action h 
appears as a new parameter. The pressure is now no 
longer uniquely fixed by dimensional analysis since in or- 
der to determine the four exponents of "f Ul a,T and h only 
three relations (from the comparison of the powers of t, 
L and M) exist. Thus the equation of state will depend 
additionally on a dimensionsless parameter x v which in- 
cludes h as well as T, a and j u . Dimensional analysis 
shows that the only possible choice isil 




Here Xt,v denotes the thermal de Broglie wave length of 
the (quasi-) particles. Ac is the Compton wave length 
in the case of massive particles. We ignore all numeri- 
cal pre-factors since they are beyond the accuracy of the 
dimensional analysis presented here. Since for h — > 0, 
i.e. Xt.v <C a, we have to regain the classical result 
(12. 4|) . quantum behavior will be seen for x v < 1, i.e. for 
a < Xt.v The borderline x v ~ 1 defines a characteristic 
energy "f u (h/a)" which turns out to be the Fermi energy 
for fermions and the temperature of Bose condensation 
for bosons, respectively. This will be discussed in detail 
in the subsequent section. 



a=X T>1 !g ( T/ mc 2 ) 




> lg (a/Xc) 



a = X T 2 

FIG. 1: The various regions (i) - (iv) of the parameter space of 
massive particles depicted in the log(T/mc 2 ) versus log(a/Ac) 
plane. Massive matter under ordinary conditions (T <C mc 2 ) 
is located in a region well below the horizontal axis. For more 
explanations see the text. 

In addition to the cross-over from the classical to the 
quantum regime, massive particles also show a cross- 
over from non-relativistic to relativistic behavior. If the 
mean particle distance is given by a, then for a 3> Xc, 
the typical particle energy -Ej p i,^/j/ approaches the non- 
relativistic expression h 2 /(ma 2 ) whereas for a<Ac the 
relativistic expression ch/a applies. 

From the mutual order of the three length scales a, Ac 
and At,i/ six regions can be distinguished 1 —. Anticipat- 
ing the result of the detailed discussions of the following 
sections these six regions exhibit very different physical 
properties : 

(i) Ac < At.2 < a : This is the region of classical be- 
havior. The only relevant length scale is the the 
inter-particle distance a. 

(ii) Ac < a < Xt,2 ■ The region of non-relativistic 
quantum behavior. Fermions are degenerated, 
bosons condense. Since Ac < a, Xt. v relativistic 
effects do not play a role. 

(iii) a < Ac < At,2 : A region of quantum behavior 
where bosons behave non-relativistically (as in (iii)) 
whereas fermions exhibit relativistic effects. The 
Pauli principle makes the difference. 

(iv) a < Xt.i < Ac : Fermions show still low temper- 
ature behavior (as in (iii)), bosons exhibit particle- 
antiparticle pair production. 

(v) At,i < a < Ac : Fermions and bosons are de- 
scribed by the same equation of state as photons. 

(vi) Xt.i < Xc < a : Same behavior as in (v). 

Bose condensation as well as particle- antiparticle pair 
production will in general change the effective inter- 
particle distance from a to At,2 and At,i, respectively. 
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Since pair production keeps some of the quantum num- 
bers constant (e.g. the total charge or the lepton or 
hadron number) in these regions a is a measure of the 
corresponding quantum number of the gas. 

For charged fermions dimensional analysis gives a fur- 
ther intrinsic length scale, the Bohr radius aeohr = 
fi 2 / (me 2 ). For simplicity we consider again a gas of elec- 
trons on a positively charged background. To have still 
an ideal gas the Coulomb interaction ~ e 2 /a has to be 
small compared to the kinetic energy. For non-relativistic 
electrons this results in e 2 /a -C h 2 /(ma 2 ), i.e., a <C 
fi 2 /(me 2 ) — aBohr- For relativistic electrons one obtains 
in the same way e 2 /a -C hc/a or e 2 /(hc) w 1/137 <C 1, 
i.e., the ideal gas condition is always fulfilled. For non- 
relativistic charged bosons outside the condensate, a has 
to be replaced by At,2 (see below) which leads to the 
condition At,2 < aBohr- The effect of Coulomb interac- 
tion will be more pronounced in the condensate since the 
kinetic energy is zero. Indeed, it has been argued that 
there is no Bose condensation in the presence of Coulomb 
interaction^. 

The different regions of the parameter space can be 
visualized in the diagram depicted in Figure ^ where we 



plot T/r 



(Ac/At.i/) 1 ' versus a/Ac in a double loga- 



rithmic representation. The lines T = mc 2 and a = Ac 
correspond to the horizontal and vertical axis, respec- 
tively. The line a = separates the classical region 
a > ^t,u from the quantum region a < Ay v . In our 
double logarithmic plot this line is straight with a slope 
equal to —v. Ordinary massive matter at temperatures 
realized on earth is restricted to the region below the 
horizontal axis. For electrons m e c 2 ~ lO 10 ^, and hence 
for room temperatures T ~ 10 2 K this region is given 
by T/mc 2 < 10 -8 . The Bohr radius as corresponds to 
as/Ac = hc/e 2 « 137. 



C. Grand Canonical Description 



appears as the only length scale of the problem, w (y) 
is so far an arbitrary but dimensionless function of fi/T. 
Following our previous dimensional arguments the pres- 
sure can now be written as^i^ 



T 



-n 



(2) 



(2.9) 



with (^) likewise so far unknown. For particles with 
vanishing chemical potential the pressure is always given 
by equation (|2.9|l with fi(0) ~ 1. Important examples 
are photons, magnons in anti-ferromagnets and acousti- 
cal phonons (all v = 1), magnons in ferromagnets (y = 2) 
and ripplons (y = 3/2). 



With (J22I) we get J = -TO(f )VX^ V and hence 
V 



x a Tr \tj t n 



(2.10) 



Using the classical equation of state l|2.4|l this requires 
dfl c i /dx = Q c i . Hence for a classical gas 



fir 



T 



fl cl {0)e^ T . 



(2.11) 



(|2.9|) and l|2.11|l can be used to calculate the chemical 
potential of the classical gas 17 



fi = Tin 



™ c; (o) 



= Tin 



a d n cl (0) 



(2.12) 



III. EQUATION OF STATE 
A. Weak quantum corrections 

Since for x v — > oo the system has to approach the clas- 
sical limit we make the following Ansatz for the pressure 



So far we assumed that the particle number N = V/a d 
and hence a was fixed. In some cases it is more appro- 
priate to consider the grand canonical description where 
instead of N the chemical potential fi is given. This 
description is in particular useful in those parameter re- 
gions where particle creation and annihilation processes 
(like the inverse /3-decay or electron-positron pair cre- 
ation) become important. 

The corresponding thermodynamic potential 
J(T, V, /i) is related to the pressure by the Gibbs- 
Duham relation 

J(T, U, M ) = VJ(T, 1, fi) = Vj(T, n) = - P V. (2.7) 

Since T and [i have the same dimension there is no intrin- 
sic length scale and hence no purely classical expression 
for the pressure. If we include fi as a further parameter 

At,, cj (£) (2.8) 



T 

P=-jiPf/b(xu) (3.13) 
a a 

with iPf/b(xu — * oo) — ► const. Here the subscript F or 
B stands for fermions or bosons. The scaling function 
iI>f/b( x v) describes the cross-over from classical to quan- 
tum behavior of particles with a dispersion relation with 
fixed value of v. For massive particles this Ansatz can be 
generalized to 

*= (3 ' 14) 

which describes in addition also the cross-over from non- 
relativistic to ultra-relativistic behavior. The scaling 
function \&f /b (^i > y) can be related to l|3.13|l in the lim- 
iting cases y — > oo and y — > 0, respectively. In the non- 
relativistic limit a 3> Ac, i.e. y — > oo, ^f/b{ x 1iU) ~ 
^F/B^xiy) = iPf/b(x2)- In the ultra-relativistic case 
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a < Ac, i-e. y — > 0, * F / B (^i,y) ~ ^f/b(&i)- In the 
following we will mainly concentrate on the discussion of 

V'F/B 

We consider first the case of weafc quantum fluctua- 
tions. Assuming that the pressure can be expanded in 
powers of the particle density, i.e. in powers of x~ d , (the 
classical result (|2.4|) is the lowest order term of this ex- 
pansion) iPf/b{%v) can be written in the form 



?Pf/b(xv) ~ 1 ± const. x v d , x v ^>\ 



(3.15) 



Here the plus and minus sign corresponds to fermions and 
bosons, respectively. Indeed, for fermions the Pauli prin- 
ciple will lead to an increase of the pressure with respect 
to the classical case. For bosons the pressure is reduced 
since quantum mechanics increases the probability for the 
double occupancy of a stated. I|2.5|l . (|3.13|) and (|3.15|) 
describe indeed the lowest order quantum correction to 
the classical equation of stated. 



B. Strong quantum limit - fermions 

In the strong quantum limit x u — > 0, we have even more 
pronounced differences between fermions and bosons. 
For fermions the Pauli-principle guarantees a non-zero 
pressure even for T — > 0. From l|2.5[) and (|3.13|) this 
requires xpp(x v ) ~ x~ v and hence 



Pf{v) 



fr v lv E P , v (a) 



1 d+v 



(3.16) 



On the r.h.s. we have introduced the Fermi energy 



E F ,u{a) 



(3.17) 



which is the kinetic energy of the fermions with momen- 
tum |p| ~ h/a. Apparently, in (|3.16|) the Fermi energy re- 
places the temperatures of the classical expression <|2.4ll . 
In particular this gives for relativistic particles^ 



Py(v = 1) 



and for non-relativistic particles 



he 



p F (u 



-,d+2 ' 



(3.18) 



(3.19) 



As long as Ac < a < At,2 (region (ii) in Figure 
1) the particles show non-relativistic quantum behav- 
ior. This situation exists e.g. for electrons in solids 
(a ~ 10 _10 m, Ac ~ 10~ 13 m) or neutrons in neutron 
stars (a ~ l(r 13 m, A c ~ lCr 16 m). 

On the other hand, for a < Ac the kinetic energy is 
larger than mc 2 and the particles behave relativistically. 
One has to expect that this 'low temperature' relativistic 
quantum behavior is seen as long as the temperature is 
small compared to the Fermi energy Ep v {a) ~ ch/a, i.e. 
for a < At,i ( corresponding to regions (iii) and (iv) of 
Figure 1). Relativistic electrons exist e.g. in white dwarfs 
(a ~ lCr 13 m). 

Ig (T/mc 2 ) 




FIG. 2: The pressure of an ideal gas of massive fermions in 
different parameter regions. The cross-over lines T = Ef,2 
and T = Ef.i separate the classical (p ~ Ta~ d ) from the 
non-relativistic quantum region (p ~ h 2 / (ma d+2 ) and the ra- 
diation (p ~ T d+1 /(hc) d ) from relativistic quantum region 
(p ~ hc/a d+1 ), respectively. Ep >v {a) denotes the Fermi en- 
ergy. The vertical line a — Ac , T < mc 2 separates the 
non-relativistic from the relativistic quantum behavior. The 
radiation region is characterized by pair creation processes 
where the effective particle distance is Xt.i < a. The sym- 
bols ♦ and it stand for the protons and electrons, respectively, 
which change their coordinates in this figure during the evo- 
lution of stars and the toy universe, respectively, considered 
in section IV. 



respectively. Both cases can be combined in the scaling 
formula 



he *f(2i, 



Pf 



Ac - 



,rf+l 



ip F (xi) 



(3.20) 



The precise form of the scaling function ^f(x, y) remains 
unknown in this approach. The cross-over between (|3.19|) 
and H3.18fl takes place at a ~ Ac, as follows from the 
general properties of ^F(x,y) discussed before. 



For very dense systems of atoms, where nuclei and 
electrons form a plasma with a <ti aeohr, electrons and 
protons may recombine to produce a neutron and a neu- 
trino (the inverse /3-decay). The charge of the resulting 
nucleus is then reduced by one. In this case the electron 
density is no longer an independent variable but deter- 
mined by the chemical potential fi of this process. If 
T <C jj. and initially a < Ac the system is in the ultra- 
relativistic region, thermal fluctuations can be neglected. 
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fl determines then unambigously a new length scale 



he 



(3.21) 



corresponding to oj(fx/T) ~ T/fi in (|2.8|> which replaces 
a in the region a < X^ < Xt,i in ea. H3.18|) . Since thermal 
fluctuations are irrelevant and a is not longer fixed, T is 
replaced by fl in all expressions. The pressure is therefore 
reduced to— 



Pf 



fi d+1 

(hc) d 



(3.22) 



We come now to a discussion of the parameter region 
of ultra-high temperatures of fermionic systems where 
particle-antiparticle pair production becomes important 
(regions (v) and (vi) of Figure 1). In this region the par- 
ticle number is not conserved and it is more adequate to 
describe the system by specifying a chemical potential, 
as we have done in subsection 111 CI For simplicity we 
assume that only fermions of the same kind (and their 
corresponding antiparticles) are generated. As discussed 
in Subsection IIIBI the Coulomb interaction can be ne- 
glected since the effective particle distance turns out to 
be much smaller than the Bohr radius as- For temper- 
atures T 3> /i the pressure is determined by l|2.9|l with 
u>(ji/T) ~ 1 and v = 1 (since we are in the relativistic 
region), i.ei 2 ^ 



T 



PF 



X 



TA 



J<d+l 



(3.23) 



region the density of particles and anti-particles, respec- 
tively, is given by 



1 / ±jl« — mc 2 

ex P 



y T,2 



T 



(3.24) 



Here we have included explicitly the term —mc 2 resulting 
from the shift of the energy scale discussed in [I] If the 
density of generated particles is small compared to the 
initial density a~ d such that X T d 2 e^ mc ^i '( T ) » a~ d we 
get after elementary manipulations for the density of the 
generated particles 



X 



2d 
T.2 



exp 



-2mc i 
T 



(3.25) 



The cross-over to the radiation region sets in if the den- 
sity (|3.25|l of generated particles is a finite fraction of 
the initial density a~ d . This gives for the cross-over line 
a co (T) (compare Figure 2) 



' mc 2 mc 2 

Xc \— eM Yd y 



(3.26) 



In the relativistic region a < Xc essentially the same 
arguments apply but now we have to replace At,2 by 
Xt,i in H3.24fl . The cross-over follows then at Xt.i < 
aexp(— mc 2 /Td) which corresponds for a < Xc essen- 
tially to a ~ Xt,i (compare Figure 2), in agreement with 
our previous conjecture. 



C. Strong quantum limit - bosons 



which agrees with the result for photons (compare 
Hence we will call (|3.23|l radiation behavior. Xta has to 
be considered as the mean distance between the fermions 
which is here smaller than a because of pair production. 
If several kinds of fermions exist each of it gives the same 
contribution l|3.23|l to the pressure. 

Finally we consider the cross-over between (|3.23fl and 
the pressure formulas valid at lower temperature. We 
begin with the region Xc < a. Lowering temperature 
we expect a cross-over to the classical result (|2.4[). To 
describe the cross-over we have to understand the contri- 
bution of pair production if we raise temperature starting 
in the classical region. It is therefore convenient to use 
a grand canonical description also for the description of 
this cross-over. To make contact with the classical re- 
sult (|2.4() . we have to express the chemical potential of 
particles and anti-particles in terms of the parameters of 
the system. Particles and anti-particles can annihilate 
into photons. Since the sum of the chemical potentials is 
conserved in this reaction 2 ^ and the chemical potential 
of photons is zero, we find that the chemical potentials 
of particle and anti-particles have opposite sign, say ±/i. 
From (|2.9|) and H2.11|) we conclude that in the classical 



For bosons the low temperature pressure cannot de- 
pend on the volume since all bosons may accumulate 
in the state of lowest energy. This requires in (|3.13() 
V'-B^y) ~ xf, from which we get for the pressure 

rp rpl + d/v / rp \d/v rp 

PbH ^ __^f_) _ (3.27) 

A T,L> 



d J v \Tc,» 



where on the r.h.s. we introduced the temperature of 
Bose condensation 



In particular, we get for relativistic particles- 5 

rp rpd-]-l 

and for non-relativistic particles 26 



(3.28) 



(3.29) 



Pb(v = 2) 



T T< d+2 )/ 2 



.d/2 



4,2 hd 



(3.30) 



6 



Again both the relativistic and the non-relativistic case 
can be combined into the scaling form 



PB 



A T.l 



■0b(zi 



(3.31) 



The cross-over between (|3.30|) and (I3.29|) occurs at Ax 1 ~ 
Ac, i-e. at T ~ mc 2 . 

Since a < Xt,v the pressure is reduced with respect 
to the classical result (|2.4p . This can be explained by 
assuming that a finite fraction of particles sits now in 
the state of zero energy and hence does not contribute 
to the pressure. Comparing (|2.4() and l|3.27|l we conclude 
that %t,v has to be considered as the mean distance of 
particles above the condensate 2 ^. 

Ig (T/mc 2 ) 




FIG. 3: The pressure of an ideal gas of massive bosons in 
different parameter regions, as discussed in the text. The line 
of Bose condensation transition T c ^{a) in the non-relativistic 
regime is continued in the relativistic regime as T Ci i(a). There 

is an additional cross-over at T ~ mc 2 from non-relativistic 
to relativistic behavior. 

Next we discuss the phase diagram of massive Bose 
particles obeying the full dispersion relation jl.ip . We 
begin with the classical region Ac < At.2 < a where the 
pressure is given by l|2.4|l . Reducing T we reach at a « 
At,2 the region where a macroscopic fraction of particles 
does not contribute to the pressure since it is in the state 
of zero energy. As argued above, the number Ne>o of 
the remaining particles of finite energy can be written as 



V 



A 



N 



T.2 



^T,2 



N 



T 



d/2 



N. 



E>Q- 



(3.32) 



Thence the number of particles in the condensate is 



A = N - N E>0 =N[1 



T 



T, 



r.2 



d/2\ 



(3.33) 



which shows that T C) 2(a) is indeed the transition tem- 
perature of Bose condensation^. If one enters the con- 
densed phase by lowering T at fixed a, the total number 
of bosons remains fixed whereas the mean distance of the 
particles above the condensate is equal to At,2 (T) , which 
increases if we further decrease temperature. 

On the other hand, if one approaches the transition 
line by decreasing a at fixed T, the mean distance At,2 
between particles above the condensate remains constant 
(since T is constant) whereas the excess particles go to 
the ground state. This behavior is seen also for a < X c 
since the effective distance of particles contributing to 
the pressure remains At,2(> Ac) (compare Fig. 3). Thus, 
Bose condensation as described by equations (|3.3()(l and 
is seen both in regions (ii) and (iii) of Figure 1. 



At sufficiently high temperatures we have to expect 
that production of particle-antiparticle pairs sets in. For 
simplicity we will again assume that only bosons of the 
same type can be created (and their anti-particles). At 
temperatures T 3> mc 2 bosons behave relativistically 
and the pressure is given by (|3.29|) . 

The discussion of the cross-over to (|2.4|l if a > Ac and 
to i|3.30[l if a < Ac, respectively, follows essentially the 
same lines as in the case of fermions, with the impor- 
tant difference that in describing the cross-over in the 
region a < At,2 the effective particle distance is now At,2 
(instead of a as in the case of fermions) . From (|3.25l) we 
therefore conclude that the cross-over occurs at T ~ mc 2 . 
Note, that the cross-over is smooth, since at T ss mc 2 , 
At,2 ~ At,i if a < Ac- 

Since in the relativistic region (iv), a < At.i < Ac, 
the effective particle density A^ 1 is smaller than a~ d one 
is inclined to conclude that here both pair creation and 
Bose condensation takes place. For a > Ac the cross-over 
is the same as for fermions. 



IV. ILLUSTRATIONS 
A. Short history of a star 

An important application of the various equations of 
state can be found in the history of stars which we 
present here in a cartoon like picture. Initially the star 
burns hydrogen to helium and other light elements and 
keeps in this way its temperature T . The ratio T/(m p c 2 ) 
is of the order 10~ 6 , where m p is the proton mass. This 
is the present situation of our sun. The pressure of a star 
can roughly be described by the classical equation (I2.4II . 
which counteracts the gravitational pressure 



PG 



m p M 

LrJV 1 „ ~ Cj 



4/3 



a 3 R 



N " 



(4.34) 



Here R and M denote the radius and the mass of the star, 
respectively, and Gjv is the gravitational constant. We 
do not present the derivation of (|4.34|) , which is straight- 
forward, but remark that pa • R 3 is the gravitational en- 
ergy of a sphere ~ GnM 2 / R, as one obtains also from 
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a dimensional argument. The ratio of the mass to the 
radius of a star in this state is controlled by the thermal 
equation of state J22J- Indeed, from the equality of (|2.4|) 
and (|4.34l) follows M/R ~ T/(Gniti p ), which increases 
linearly with temperature. 

Once the nuclear fuel is exhausted, temperature can- 
not be kept and the star collapses under the graviational 
force. As a result the density increases and the star 
matter crosses over to the region (ii) of Figure 1, where 
the pressure is given by (|3.19(1 . Since the pressure is 
now inversely proportional to the mass, mainly the elec- 
trons contribute to it. Once under further contraction a 
has reached Ac (corresponding to densities of the order 
~ 10 7 gcm~ 3 ) (|3.18|) has to be applied for the pressure. 
The star has reached the state of a white dwarf (compare 
Figure 2). 

If the density is increased further ( i.e. a < A^), neu- 
trons are generated via the inverse /3-decay. The pressure 
is determined then by (|3.22|) until most of the nucleons 
became neutrons. For the pressure in neutron stars for- 
mulas (|3.19|) for a > Xc and (|3.18|) for a < Xc again 
apply, but the mass m„ ss m p is now that of the neu- 
trons (compare Figure 2). This is the fate of stars with 
a mass between 15 and 30 times the mass of the sun. 

In the ultra-relativistic regions of white dwarfs and 
neutrons stars the equality between (|4.34(l and (|3.18|l re- 
sults in a relation for the mass 

This so-called Chandrasekhar mass Mq is the maximum 
mass of white dwarfs and neutron stars. Its value is de- 
termined by the thermal equation of state in region (hi) 
and is of the order of twice the mass of the sun. 



B. A toy universe 

As a second illustration let us consider a contracting 
toy universe . For simplicity, this universe is supposed to 
have no gravity, i.e. there are no stars, clouds of cosmic 
dust etc., just photons and diluted massive particles. To 
make things even simpler the only particles assumed to 
exist are hydrogen atoms ( a "non-metallic" universe in 
the language of astronomers). Initially the photon spec- 
trum is assumed to correspond to black body radiation 
of about 3K. Since for massless particles At,i is the only 
existing length scale, the typical distance of photons is of 
the order Xt.i ~ 10 -3 m. 

The number of H-atoms per m 3 is assumed to be of 
order 1. Hydrogen atoms are bosons and hence could 
be represented in Figure |3J The horizontal and vertical 
coordinates of the hydrogen gas in Figure |21 are given by 
o-h/^ch ~ 10 16 and T/mc 2 « 10~ 16 , respectively, i.e. 
they are outside of the margins of this figure. Xq.h ~ 
2.1 ■ 10~ 16 m is the Compton wave length of hydrogen. 



Next we switch on the cosmic contraction, i.e. we go 
backwards in time. It is convenient to introduce the con- 
traction factor r(t) = a(t)/a(tf) < 1, where a(t) de- 
notes the distance at an earlier time t < tf when the 
final distance was a(tf) = a. Since contraction of space 
reduces the wavelength of photons such that At,i(^) = 
*T,i(*/)K*) = hc/T^t) with T 7 (i) = T 7 (i/)/r(t), i.e. 
the photon gas heats up. 

Without coupling to the photons, the hydrogen atoms 
form a classical non-relativistic gas (region (i) of Fig- 
ure 1) which contracts adiabatically. Its temperature 
ThU) = TH(tf)/r 2 (t) rises proportional to r~ 2 (t), pro- 
vided there is enough interaction between the hydro- 
gen atoms to reach thermal equilibrium. The horizon- 
tal and vertical coordinates of the hydrogen gas in Fig- 
ure 3 thence change according to an{t) / Xc, P ~ r(t) and 
T(t)/mc 2 ~ r~ 2 (t), respectively. 

When the temperature of the photons reaches about 
3000 K (r(t) ~ 10~ 3 ), the hydrogen atoms have a dis- 
tance of about a nit) ~ 10~ 3 m and are now par- 
tially ionized, i.e. we have free protons and electrons, 
both fcrmions, which appear in Figure El Since both 
particles are charged, they will emit and adsorb pho- 
tons. From now on matter and radiation have approx- 
imately the same temperature T(t) ~ ?>r~ 1 {t)K. The 
initial coordinates are a p (t)/Xc P ~ 10 13 , T/m p c 2 = 
&C, v /*t,u) v ~ 10- 10 ( i.e. a p (t)/X T .2 p ~ 10 8 ) for pro- 
tons and a e (t)/X c ,e ~ 10 10 , T/m e c 2 = (X c ,d*T,v) v ~ 
10~ 7 ( i.e. a e (i)/Ar,2 e ~ 10 6 ) for electrons, respectively. 
Thus their initial positions are in region (i) of Figure 
1. As long as a T /Xc ~ 10 4 -C a(t)/Xc, electrons and 
protons form initially an non-interacting non-relativistic 
classical plasma (compare Figure 2). Since a(t) ~ r(T) 
and T(t) ~ r -1 (i), the slope of the path under further 
contraction in the double logarithmic plot of Figure 2 is 
— 1 and remains unchanged until we reach the relativistic 
region T > mc 2 . 

When the contraction parameter reaches r(i) ~ 10 -10 , 
i.e temperatures reach about 3 ■ 10 10 i<r, the mean particle 
distance is of the order of the Bohr radius. 

When temperatures reach the rest energy of the elec- 
trons, i.e. crossing the bisecting line T = m e c 2 of 
Figure 2, the electrons start to behave relativistically. 
For T m e c 2 electron-positron pairs are formed and 
leptons behave ultra-relativistically. For temperatures 
T > lO 15 ^ and a p ~ 10 _15 m hadrons decay into quarks 
and anti-quarks, which also behave ultra-relativistically. 
Also their mean distance is now determined by pair pro- 
duction. In addition gluons, i.e. the bosons which carry 
the interaction between quarks, have to be considered. 
Eventually the coordinates of all fundamental fermions 
and bosons reach the line a ~ At,i- 

The pressure in all cases is that or ultra-relativistic 
fermions (|3.23|l or bosons l|3.29|l . respectively. The total 
pressure is the sum of these contributions, i.e. it is pro- 
portional to the number of different fermions and bosons. 
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V. ADDITIONAL THERMODYNAMIC 
RELATIONS 



We can also use the results of the previous Chapters to 
determine the scaling behavior of the free energy f(T, v) 
per particle which can be written as 

f( T , v ) = - Ti Pf/b(xv) ■ (5-36) 
Differentiation with respect to v — a d gives 
df _ Tip' p/B [x u )x u _ 



dv 



vd 



(5.37) 



Thus iPf/b(xv) = 0-/d)x u cp' F , B (x u ). The entropy per 
particle s follows from 



and hence we get for the energy per particle u 

v 1 v 



(5.38) 



(5.39) 



which is the well known relation between the energy u 
and pv both for relativistic and non-relativistic particle in 
any dimension^. For the chemical potential fj, we obtain 
from the Gibbs-Duham relation 

H = f+pv = T \^<p F/B (x v )x u - ip F/B (x, y )^j . (5.40) 

With our Ansatze used above we find at low temperatures 
s = for fermions and fi = for bosons. 

These expressions can be used to calculate additional 
thermodynamic quantities. The specific heat follows from 
(|5~57|) as 

CT = Wf = = vW^*^ ■ (5 ' 41) 

For fermions our Ansatz ipp(x v ) ~ x~ v ', i.e. / ~ Ep, 
gives ct — instead of the correct result ct ~ T/Ep ~ 
x". This is not surprising since we determined vp ~ / 
from the condition that the free energy should go to a 
constant (the Fermi energy) for T — > 0. 



To get the correct specific heat we have to add finite 
temperature corrections to u. Indeed, the low tempera- 
ture corrections to energy of the Fermi gas result from 
the difference A of the Fermi distribution at non-zero 
and zero temperatures (the latter is the step function). 
These corrections are of the formS 



dE g(E) E A 



E 



T 



(5.42) 



where A(x) = (e x + 1) _1 — 9(— x) , and g{E) is the 
density of states which is smooth at E = \i. With 
A(x) = — A(— x) it is clear that the low temperature 
expansion of u includes only even powers in T. Since Ep 
is the only energy scale, the low temperature expansion 
of the energy of the ideal Fermi gas is therefore of the 
form 



T 2 

M1 + 33T 



rp4 



(5.43) 



which is the Sommerfeld expansion. This gives for the 
specific heat ct ~ T/Ep2&. On the other hand, cal- 
culating the compressibility we get with ^}p(x v ) ~ x~ v 
immediately the correct result 



K = — — 



\_dV_ 
V~dp 



For bosons we obtain with ipg 



(5.44) 



x d the correct result 



for the specific heat^ c T ~ ^(Tx d ) - {T/T c ) d,u 



(«/A, 



and the compressibility diverges. 



To conclude we have shown that the main properties 
of ideal quantum gases can be obtained from simple di- 
mensional arguments and the Pauli principle, without 
resorting to statistical mechanics. 
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